The Mott- Anderson transition in the disordered one-dimensional Hubbard model 
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We use the density matrix renormalization group to study the quantum transitions that occur in 
the half-filled one-dimensional fermionic Hubbard model with onsite potential disorder. We find a 
transition from the gapped Mott phase with algebraic spin correlations to a gapless spin-disordered 
phase beyond a critical strength of the disorder A c ~ U/2. Both the transitions in the charge and 
spin sectors are shown to be coincident. We also establish the finite-size corrections to the charge 
gap and the spin-spin correlation length in the presence of disorder and using a finite-size-scaling 
analysis we obtain the zero temperature phase diagram of the various quantum phase transitions 
that occur in the disorder-interaction plane. 



71.27+a, 71.30.+h 

Electronic systems can undergo quantum phase tran- 
sitions from metallic to insulating behavior as a function 
of either the interaction strength or the degree of disor- 
der or both [Q . A clean system at certain commensurate 
fillings may develop a gap in the energy spectrum as the 
strength of the repulsion is increased and turn into a Mott 
insulator |Q. On the other hand, a system of noninter- 
acting electrons can undergo a transition from metal to 
insulator as the degree of randomness is increased and 
turn into an Anderson localized insulator J3j. The inter- 
play of electron-electron interaction and disorder raises 
interesting possibilities of a new type of transition, dis- 
tinct from the clean correlation induced Mott transition 
or the disorder induced Anderson transition. 

The repulsive Hubbard model in one dimension (ID) 
is probably the simplest model which shows a Mott tran- 
sition at half-filling for arbitrary values of the interaction 
strength U > B. One of the most attractive features of 
this model is that the Mott transition is unaccompanied 
by a spin-density-wave gap as e.g. happens in the 2D 
Hubbard model due to the existence of magnetic long- 
range order (MLRO) ||. Instead the ID model shows 
algebraically decaying spin correlations as the maximal 
remnant of MLRO in ID. The properties of the clean ID 
Hubbard model are well established: The dependence 
of the Mott gap on the interaction parameter [Q, the 
asymptotic behavior and the critical exponents of vari- 
ous correlation functions [|| have all been computed us- 
ing Bethe Ansatz and the finite-size-scaling approach of 
conformal quantum field theory [Q. All this makes this 
model particularly attractive to study the effect of disor- 
der on the Mott state. 

Recent studies of the half-filled disordered ID Hub- 
bard model have proceeded numerically using the Quan- 
tum Monte Carlo (QMC) method || and analytically us- 
ing bosonization and the renormalization group method 
M. The QMC results give a very accurate description 



of the finite temperature properties of the system. How- 
ever, the low temperature properties can only be inferred 
by extrapolating the finite temperature data and going 
to larger system sizes. On the other hand, bosoniza- 
tion methods using perturbative renormalization group 
techniques only give an indication of the various plau- 
sible fixed points in parameter space. Computing cor- 
relation functions at the strong coupling fixed points is 
then impossible as the coupling constants are driven away 
from their weak coupling values. Nevertheless, the re- 
sults of both approaches indicate that at half-filling a 
finite amount of potential disorder is needed to cause a 
transition from Mott (gapped) to Anderson (gapless) in- 
sulating behavior. This is in qualitative agreement with 
arguments put forward by Ma JToJ . 

In the present Letter, we have studied the ground 
state properties of the disordered ID Hubbard model 
at zero temperature with the help of the density-matrix 
renormalization group (DMRG) |nj. This method has 
been previously shown to be highly successful for ID 
quantum systems fl^ , |l3]| and may be seen as a numer- 
ical variational- wave- function approach 14 1 . After in- 
troducing the parameters of our DMRG, we show that 
in the clean case, the previously known results can be 
reproduced numerically reliably. We then consider fi- 
nite disorder and show that the charge gap G c remains 
open for small disorder up to a critical disorder strength 
A c U/2. The transition in the spin sector is seen by 
studying the behavior of the spin-spin-correlation func- 
tion (S~(r)S + (0)). For A < A c , the power-law remnant 
of the MLRO persists, whereas for A > A c , the spin-spin 
correlation indicates the emergence of a spin-disordered 
phase. 

The Hubbard Hamiltonian with additional potential 
disorder on a chain of L sites is given as 
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H = -t {c f x+lcr c X o- + h.c.) + U^2 n x-\n xl + ^2 

x=\ x=l x=l 

(1) 

where —t is the hopping amplitude between nearest- 
neighbor, c xa .(c xa ) the fermion creation (annihilation) 
operator at site x with spin er, n xa = c x<J c X(J the number 
operator, n x — n x | 4- n x | and U is the onsite repulsive 
energy. The onsite chemical potential p x is a random 
number which we take to be uniformly distributed be- 
tween ±A such that A = corresponds to the clean 
case. We work at half-filling, e.g. N — (J2 x =i n x) = L 
and the energy scale is set by choosing t = 1. 

We follow the standard open chain DMRG algorithm 
of White jnj. In each iteration, we diagonalize the 
Hamiltonian matrix of a super-block denoted by B L l * 
• B r L , 2 _ 1 of L sites and obtain the energy and the wave 
function |V>ol) of the ground state fL3f . Here • repre- 
sents a single site. Using IV'ol) as the target state Q we 
compute the reduced density matrix p l of the left sub- 
block flL = S£/ 2 -i* °f s i ze L/2. Diagonalizing p l we 
obtain its eigenstates and the eigenvalues. These eigen- 
states with the highest eigenvalues are the most probable 
states of the left sub-block when the super-block is in the 
state \iPol) and so can be used for truncation. Keeping 
M eigenstates corresponding to the largest M eigenvalues 
of p l as the new basis of the left sub-block, we transform 
the Hamiltonian and all further operators into this new 
basis. Since every single • has 4 states, we thus truncate 
the original 4M states of B l L , % to only M states. Usually 
we have used M = 128. Because the disorder destroys 
translational symmetry, the left and right sub-blocks are 
non-identical. Hence a similar procedure needs to be fol- 
lowed for the right sub-block B r L ^ 2 . The above steps are 

then repeated for the new super-block B L , 2 • • B r L j 2 with 
L + 2 sites and thus the system increases by two sites at 
each iteration. 

Let Eq(L) denote the ground state energy of the Hub- 
bard chain of length L with iVj = N± = L/2. In order to 
compute the charge gap G C L for a system of finite length 
L, we repeat the DMRG steps with Nf = L/2 + 1 and 
N± = L/2 and also with 7V T = L/2 - 1 and JVj. = L/2. 
We denote these two ground state energies by E\{L) and 
E—i(L), respectively. The charge gap is then defined as 
the discontinuity of the chemical potential at half-filling 
@, i.e., G C L = Ei(L) + E_ X (L) - 2E (L). In the pres- 
ence of disorder G C L is computed for at least 10 differ- 
ent disorder realizations and then averaged over all such 
realizations. We also compute the spin-spin-correlation 
function T s L (r) = (ip OL \S-(r)S + (0)\ipo L ) and the sec- 
ond moment of the staggered antiferromagnetic (AFM) 
correlation function (££) 2 = ^2 r 2 (— l)T£(r). Here, 
S + (r) = CrfCj-i. In the presence of disorder, the correla- 
tion functions are first averaged over the disorder realiza- 



tions and then the correlation length of the averaged 
staggered correlation function is found. 

The clean case A = 0: It has been shown in Ref. 0] 
that the half-filled repulsive ID Hubbard model exhibits 
a charge gap for all non-zero values of the interaction 
strength. We find that the functional dependence of the 
thermodynamic value of the charge gap G^ on the in- 
teraction strength U as obtained by DMRG (Fig. [I]) is in 
excellent agreement with the exact solution computed in 
Ref. @, 

G c = U - 4 + 8 r dw— Jl( - W ) (2) 

with Ji(w) a Bessel function. We note that the DMRG 
algorithm using open boundary conditions gives consis- 
tent results for the extrapolated charge gap with that 
of Eq. p| which has been derived using periodic boundary 
conditions. Nevertheless, we now employ a finite-size- 
scaling (FSS) analysis of the charge gap G C L in order to 
remove any finite-size effects that could arise in the ex- 
trapolation. In Fig. [I] (inset) we show that the leading 
order finite-size corrections to G^, fall off as 1/L, i.e. 
G C L {U) = G^(U) + g{U)/L. Having determined the ex- 
plicit scale dependence of G C L , a plot of LG C L versus U as 
in Fig. U shows curves for different L coalescing as the 
charge gap vanishes. This finite-size behavior allows a 
numerically accurate determination of the critical value 
of the interaction strength and we find in accordance with 
the result of Ref. H that the Mott transition occurs at 
U = 0. 

The disordered case A ^ 0: Quantum Monte Carlo 
H and bosonization || studies have predicted the ex- 
istence of a critical disorder A c > beyond which the 
Mott gap vanishes. In the limit of large U, this can be 
motivated p5[ by potential energy considerations: any 
rearrangement of the one particle per site configuration 
in the half-filled Hubbard model would necessarily cost 
an energy U due to double occupancy and would gain in 
the local site potential energy a maximum of 2A. Hence, 
in order for it to be feasible for the electrons to take 
advantage of the random site energies we should have 
A c > U/2 in the case of bounded disorder. In Fig. |^ we 
show G^, as a function of the disorder strength A for 
U = 2. We see that the system undergoes a transition 
from a gapped Mott insulator phase to a gapless phase 
for A c w 1 = U/2. Fig. [| (inset) shows the finite-size 
corrections to G^ in the presence of two representative 
weak (A = 0.1) and strong (A = 1.0) disorder strengths 
for U = 2. We note that the corrections continue to fall 
off as 1/L even in the presence of disorder. The FSS plot 
of LG C L versus A (Fig. ||) shows that curves for different 
L coalesce at A c = 1.0 ± 0.1 again indicating the tran- 
sition from a gapped to a gapless phase as the strength 
of the disorder is increased. We have also done the same 
analysis for U = 3 and 5 where a well developed Mott 
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gap exists and where the DMRG is more stable |Tl]] . Our 
results indicate that the transition into the gapless phase 
takes place at A c w U/2. 

Spin- spin- correlation function: For A = 0, the half- 
filled ID Hubbard model shows algebraically decaying 
antifcrromagnetic (AFM) correlations as computed by 
the methods of conformal field theory and bosonization 



r-(r)^ V os( " r + (M +o(iA- 2 ), 



(3) 



When the spin-spin-correlation function decays as a 
power-law, as in Eq. 0, the AFM correlation length £^ 
diverges in the thermodynamic limit. However, in the 
absence of long-ranged correlations, £^ remains finite. 
Fig. H (inset) shows that for A = 0.1 the finite-size cor- 
rection to continues to fall off as 1/L such that 
1/££(A) = 1/£^(A) + C(A)/L and we can also apply the 
previous FSS analysis to the correlation length. In Fig. 
I we plot L/£ S L as a function of A. The data for different 
L coalesce until the disorder A s « 0.9 ± 0.1. Thus for 
U = 2 the staggered T s (r) continues to fall off with a 
power-law up to this critical disorder A s . For larger A 
beyond A s , the values of L/^ S L do not coalesce any more. 
This indicates the transition from power-law correlations 
into a short-ranged spin-disordered phase with a finite 
correlation length. Again, the same analysis of the spin- 
spin-correlation function for U = 3 and 5 confirms that 
A, « U/2. 

Thus the two transitions seen above viz. (I) gapped 
Mott insulator phase to a gapless phase at A c and 
(II) long-ranged AFM phase to a short-ranged spin- 
disordered phase at A s , are coincident, and A c = A s 
within the numerical accuracy. Beyond the critical dis- 
order A c w U/2 the gapped Mott insulator with power- 
law AFM correlations goes over to a gapless short-ranged 
spin-disordered phase. This is the main result of our work 
and allows us to show in Fig. ^ the phase diagram of the 
zero temperature quantum transitions that occur in the 
ID disordered half-filled Hubbard model. 

In summary, the DMRG allows us to analyze the var- 
ious quantum transitions that occur in the half-filled ID 
Hubbard model in the presence of onsite disorder. For 
A = we find that the finite-size corrections to the ther- 
modynamic value of the charge gap G^ scale as 1/L 
(Fig. |l| (inset)). The functional dependence of G c ^ ) (t/) on 
the interaction strength (Fig. ^) is in excellent agreement 
with the exact solution computed in Ref. H in support 
of our DMRG approach. We have further shown that 
for A ^ 0, the finite size corrections to G^(A) continue 
to fall off as 1/L. For small A < U/2, the Mott gap 
is shown to survive. Curves of LG C L versus A for differ- 
ent L (Fig. |4|) come together at A c (w 1.0) and coalesce 
after that indicating a transition from the gapped Mott 
insulator phase to a gapless phase. A finite-size anal- 
ysis of the staggered AFM correlation length shows 



that the values of L/^ S L for various disorders A (Fig. |5|) 
also coalesce at A s « 1.0 ± 0.1. These transitions in the 
charge and the spin sector are shown to be coincident 
A c = A s . Thus in the gapped Mott phase T s (r) decays 
algebraically and is short ranged in the disordered gap- 
less phase. We obtain the (U, A) phase diagram (Fig. ||) 
showing the phase boundary which separates the gapped 
Mott insulator phase with algebraic spin correlations and 
the spin-disordered gapless phase. 



ACKNOWLEDGMENTS 

We sincerely thank S. Ramashesha for useful discus- 
sions on the DMRG technique and the HRZ (TUCZ) and 
SERC (IISc) for the computing facilities. This work has 
been supported by the Deutsche Forschungsgemeinschaft 
(Sfb 393). 



For a review, see, P. A. Lee and T. V. Ramakrishnan, 
Rev. Mod. Phys. 57, 287 (1985); D. Belitz and T. R. 
Kirkpatrick, Rev. Mod. Phys. 66, 261 (1994). 
[2] N. F. Mott, Metal-Insulator Transitions, Taylor & Fran- 
cis, (London 1990). 
[3] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 
[4] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 
(1968). 

[5] S. R. White, D. J. Scalapino, R. L. Sugar, E. Y. Loh, J. 
E. Gubernatis, and R. T. Scalettar, Phys. Rev. B 40, 506 
(1989). 

[6] F. Woynarovich, J. Phys. A 22, 4243 (1989); H. Frahm 

and V. E. Korepin, Phys. Rev. B 42, 10553 (1990); H. J. 

Schulz, Phys. Rev. Lett. 64, 2831 (1990); 
[7] A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, 

Nucl. Phys. B241, 333 (1984); J. L. Cardy, Phys. Rev. 

Lett. 52, 1575 (1984). 
[8] A. W. Sandvik, D. J. Scalapino and P. Henelius, Phys. 

Rev. B 50, 10474 (1994). 
[91 S. Fujimoto and N. Kawakami, preprint (1996, cond 



[10 

[11 
[12 

[is; 
[14; 
[15; 



mat/9609011) 



M. Ma, Phys. Rev. B 26, 5097 (1982). 
S. R. White, Phys. Rev. Lett. 69, 2863 (1993). 
For a brief review, see, G. A. Gehring, R. J. B ursill, and 
T. Xiang, preprint (1996, |cond-mat/9608127| ). 
R. V. Pai, R. Pandit, H. R. Krishnamurthy, and S. Ra- 
mashesha, Phys. Rev. Lett. 76, 2937 (1996). 
S. Ostlund and S. Rommer, Phys. Rev. Lett. 75, 3537 
(1995). 

R. Shankar, Int. J. Mod. Phys. B 4, 2371 (1990). 



3 




12 



10 



FIG. 1. Filled squares (■) show the charge gap G^ as a 
function of the interaction strength U. The exact result is 
shown by the solid line. Inset: G C L vs 1/L for U = 0.0, 2.5, 3.0 
and 3.5. The value of the intercept gives G^. 
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FIG. 4. LG C L as a function of the disorder A for (7 = 2 
showing the coalescence of the curves for different L at 
A = 1.0. This indicates the transition from a gapped Mott 
insulator phase to a gapless phase. 
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FIG. 2. LG C L as a function of the interaction strength U 
showing the coalescence of curves for different L at U = 0. 
This indicates that the critical value of the interaction 
strength at which the Mott transition occurs is U = 0. 
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FIG. 3. Filled squares (■) show the charge gap G^ as a 
function of the disorder A for (7 = 2. Inset: G C L vs 1/L for 
A = 0.1 and 1.0. The value of the intercept gives G^. 




FIG. 5. L/££ as a function of the disorder A for (7 = 2 
showing the coalescence of the curves for different L at 
A ~ 0.9 — 1.0. This indicates the transition from a gapped 
Mott insulator phase to a spin-disordered gapless phase. In- 
set: vs 1/L for A = 0.1. 




FIG. 6. The phase diagram of the disordered ID Hub- 
bard Hamiltonian (|l]) at half-filling showing the Mott insu- 
lator with algebraic spin correlations (I) and the gapless spin 
disordered phase (II). The phase boundary has been drawn 
through the computed points (filled squares (■) with error 
bars) and is close to A = U/2. 
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